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Abstract 

We study the adsorption-desorption phase transition of directed branched poly- 
mer in d+ 1 dimensions in contact with a line by mapping it to a d dimensional hard 
core lattice gas at negative activity. We solve the model exactly in 1 + 1 dimensions, 
and calculate the crossover exponent related to fraction of monomers adsorbed at 
the critical point of surface transition, and we also determine the density profile of 
the polymer in different phases. We also obtain the value of crossover exponent in 
2 + 1 dimensions and give the scaling function of the sticking fraction for 1 + 1 and 
2 + 1 dimensional directed branched polymer. 

Linear and branched polymers, near an attractive surface, undergo an adsorption- 
desorption transition, which has important applications in areas ranging from technology 
such as in lubrication, adhesion, surface protection to biology [1, 2]. For example, adsorbed 
polymers are used for surface-modification of medical implants [3]. There have been 
several theoretical studies of the behavior of a polymer near a surface[2, 4, 5, 6, 7, 8]. 
Especially, effect of surface for idealized polymer (with no self-exclusion), modeled by 
random walks has been studied extensively. There are many exact results known for 
Gaussian random walks in presence of a surface [9, 10]. In comparison, linear polymer 
with self exclusion and branched polymer are less well studied. For a self-avoiding walk 
(SAW) in the vicinity of a surface the exact critical exponents are known from conformal 
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field theory [11]. Directed polymer chain adsorption, modeled by a directed SAW is one 
of the few solvable models of surface effects in 2 and 3 dimensions [6]. For directed 
walks, self exclusion is automatic, and nontrivial effects of excluded volume interaction 
are not seen. For branched polymers (modeled by lattice animals), a relation between 
the exponent characterizing the number of animals, in presence of surface, and in the 
bulk is known from a simple argument given by De'Bell et. al [12]. In this paper we 
solve directed branched polymer (DBP) in 2 dimensions and 3 dimensions in presence of 
a Id line exactly. Introducing a preferred direction makes the system analytically more 
tractable. Similar results for 2 dimensions have also been obtained by Rensburg et. al. 
[13] independently. 

The enumeration of directed site animals in d + 1 dimensions is related to hard-core 
lattice gas (HCLG) at negative activity with repulsive interactions in d dimensions and 
the Yang-Lee edge problem in d dimensions [14, 15, 16, 17]. In this paper we give the 
mapping of a d + 1 dimensional directed branched polymer in presence of a line to a d 
dimensional HCLG with repulsive interactions. 

The plan of the paper is as follows. In section 1 we will define the model of directed 
branched polymer and the quantities of interest. Using the above mentioned correspon- 
dence we give the mapping of a d + 1 dimensional directed branched polymer in presence 
of a line to a d dimensional HCLG with repulsive interactions in Section 2. In Section 3, 
for 1 + 1 d DBP in presence of a Id penetrable line we solve the model exactly. Section 
4 deals with the DBP in 1 + 1 d in presence of a impenetrable line. For DBP in 1 + 1 
dimensions we show that the behavior at the transition point for penetrable and impene- 
trable wall is the same, and not just the crossover exponent but even the density profile is 
the same. This implies that for 1 + 1 dimensions, for impenetrable surface, at the phase 
transition point the decrease in entropy is exactly compensated by the increase in internal 
energy. This seems to be a special property of polymers in 2 dimensions. Even for linear 
polymers the exponent for both cases is the same and hence it is believed that for a linear 
polymer in 2 dimensions in presence of a impenetrable surface the phase transition point 
correspond to point where surface effects vanish completely and system behaves like bulk 
[8]. Here we are able to show it explicitly for directed branched polymers. In Section 5 
using Baxter's solution of hard hexagon gas we study 2 + 1 dimensional DBP in presence 
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of line and calculate the crossover exponent and sticking fraction for the directed branched 
polymer exactly. The scaling function of sticking fraction is a function of two intensive 
thermodynamic variables. We have derived its exact form in 1 + 1 dimension and 2 + 1 
dimension. There are very few such exact nontrivial scaling functions of more than one 
thermodynamic variable known [18]. We also get the large w expansion of the sticking 
fraction as a power series in 1/y/w. 

1 The Model 

A directed branched polymer or a directed animal on a lattice, rooted at the origin is a 
connected cluster such that any site of the animal can be reached from the root by a walk 
which never goes opposite to the preferred direction. For example, on a square lattice 
drawn tilted at 45% in Fig. 1, a directed site animal or a directed branched polymer A 
rooted at the origin is a set of occupied sites including origin, such that for each occupied 
site (x, t) other than the origin, at least one of the two sites (x — l,t — 1) and (x + l,t — l) 
is also occupied. The number of sites in A will be denoted by s = \ A\. We define n(x\A) 
as the number of sites of A having the transverse coordinate x. We study the DBP in 
presence of Id line parallel to the preferred direction. This is positioned along the main 
diagonal of the lattice (Fig. 1). We will consider only polymers rooted at the surface in 
this paper. 

We assign a fugacity y to all allowed sites of the cluster. Further, if we associate an 
additional energy — E with each site on the surface, each site on surface will have an 
additional weight and the fugacity of sites about the diagonal, denoted by yo is equal to 
wy where 

w = exp(E/kT) (1) 

Hence w > 1 would correspond to an attractive surface. 

We define A(w,y), the grand partition function of the polymer as 

oo 

A(w, y) = £ yWw n ° = ]T A s (w)y s (2) 

A s=l 

where no = n(0|*4.) and A s (w) is the partition function of the polymer made of exactly s 
monomers. 
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Figure 1: A Directed Branched Polymer of size 50, rooted on the surface. 

For w — 1, we get the statistics of equally weighted animals and A s (l) is the number 
of distinct directed animals having s sites with given boundary conditions. For large 
s, A s (w) varies as X s s e , where 9 is known as entropic critical exponent. Similarly, the 
transverse size of the polymer for large s scales as s u , where v is the exponent which 
defines the transverse length scale of the polymer. These exponents take different values 
in desorbed , adsorbed and in the critical regions. We will use subscripts de, c and ad 
to represent critical exponents and other quantities in desorbed , critical and adsorbed 
phases of the polymer. 

And the free energy per monomer of the polymer in thermodynamic limit is given by 

F(T) = hm — logA s (w) = k B Tlog(y 00 (w)) (3) 

where Uoo{w) is the value of fugacity at which A(w, y) has a singularity for a given value 
of w. 

Let 0(x, s) be the value of n(x.,A) averaged over all configurations A of size s. We 
define a generating function ty(x.;w,y) as 

*(x; w,y) = J2 n(x\A)w n °y^ = £ 0(x, s)A.(w) y s (4) 

A s 



4 



There is a critical value w c of wall activity such that for w > w c , 0(0, s) is proportional 
to s for large s and the transverse size is finite {v a( i = 0). This is the adsorbed phase, 
in which monomers tend to stick to the surface, w < w c corresponds to the desorbed 
phase of the polymer in which only a finite number of monomers stick to the surface. At 
w = w c , the critical point of the surface transition, the number of adsorbed monomers as 
function of polymer size in large s limit have a behavior given by 

<j) c (0,s)~s a ; (5) 

where a is known as the crossover exponent of the surface transition. 

In the s — > oo limit, the fraction of monomers adsorbed is like order parameter of 
the surface phase transition. In the constant fugacity ensemble A(w,y) is the partition 
function with fixed w and y and hence the average polymer size would be given by 

/ / n X Esy s w no d\nA(w,y) 

(sly.w)) — — = (6) 

Similarly, the average number of monomers at the surface would be 

Ey s w no d\nw v ; 

The sticking fraction defined as the fraction of polymer segments at the surface, rep- 
resented by C st (w,y), would be given by 

(n {y,w)) 

C st {w,y) = — 8 

(s(y,w)) 

In the infinite polymer limit, if we represent the value of fugacity at which (s(y,w)) 
diverges by yoo(w) for a given w, then the sticking fraction is only a function of the wall 
activity w and is given by 

dln yoo {w) 

Cst{w) = —a^r (9) 

This is the order parameter of the surface phase transition and is zero for w < w c , 
where w c is the surface phase transition point. 

In general, in the large polymer limit, near critical value of w, as w — * w+, C st (w,y) 
is expected to have scaling form 
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C st (w,y) = e 1 - a h((w-w c )e~ a ) (10) 

where e = 1 — y/y^w). The scaling function h(u) where u = (w — w c )e" a , is a function of 
w and y which are both intensive thermodynamic variables. As u — > oo, ~ u( 1-a )/ a . 

2 General Results 

The directed site animal enumeration (DSAE) problem in d + 1-dimensions is related to 
time development of thermal relaxation of a hard core lattice gas (HCLG) with nearest 
neighbor exclusion on d dimensional lattice [15]. In [19], we have shown that this corre- 
spondence relates the density at a site i in steady state to sum of weights of all animals 
rooted at i, the grand partition function of the animal. Also, the average number of sites 
at a given transverse distance x from the origin for a d + 1 dimensional directed animal 
is related to the density-density correlation function of the lattice gas in d dimensions. 

Specifically, if on a d+1 dimensional body-centered hyper-cubic lattice we define weight 
of an animal A as the product of weights of all occupied sites, with weight corresponding 
to a site with x coordinate i being y^ then the DSAE problem on this d+1 dimensional 
lattice gets related to time development of HCLG with nearest neighbor exclusion on 
a d dimensional body-centered hyper-cubic lattice with the rates which satisfy detailed 
balance condition corresponding to the Hamiltonian 

H = +oo n i n -j - XX m ^) n * ( n ) 

<ij> i 

where Zi = — yi/(l +yi) and the animal number generating function is just the negative of 
density of HCLG with change of variables from z to y. Here we have used the convention 
that if J2<ij> n i n j — then the corresponding term in the Hamiltonian is zero. The 
configurations with any pair of occupied nearest neighbor have infinite energy and do not 
contribute to the partition function. 

The partition function is linear in all z[s. The linearity of the partition function in 
z\s implies that in case when the activity about x = is different from that in the rest of 
the sample, i.e if we let the activity about x = be z and activity in rest of space be z, 
then the partition function of the HCLG can be written as 
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Z(z ,z)=A(z)+z B(z) (12) 

where A(z) and B(z) are polynomials in z. If p represents the density of HCLG when the 
activity about each site is the same, then the density of HCLG about the origin in the 
present case po(z , z) can be written in terms of p as 

Po(zo,z) = -^j- r (13) 

pz + z(l - p) 

Same observation has been made by Cardy in [20]. Correspondingly, since A(w,y) is 
just the negative of po(zo,z) with zo = —wy/(l + wy) and z = —y/(l + y), we can express 
A(w,y) in terms of A(l,y) and this is given by 

1 ,y) (l + wy)+A(l,y)(l-w) 1 ; 

Moreover, the density- density correlation function of HCLG G(x; w, z) with tu^l can 
be expressed in terms of density density correlation function when w — 1. We find that 
the density density correlation function is related to \I/(x; y) on a hyper-cubic lattice 
as follows 

tf(x;iu,w) = — G fx: w, z = — ^- ) (15) 



from this we get 

V(x;w,y) _ w(l + y)[l + wy-(l-A(l,y))(l-w)} 



(16) 



*(x;l,y) [l + wy + A(l,y)(l-w)} 2 

Since p is the density of the HCLG, then as discussed in [19], for x = 0, the density 
density correlation of HCLG is always equal to p(p — 1) for any d dimensional case and 
hence \l/(0; l,y) can be completely expressed in terms of A(l,y). Hence we get, 

mv>,v) = * 1+mi y i \t A{1 -$! (17) 

ll + w y + A(l, y )(l-w)]' 
Eq. (14-17) hold for all dimensions. Hence, in presence of Id surface, a DBP in d + 1 
dimensions rooted on the surface can be studied using the mapping to HCLG. Moreover 
the generating functions A(w,y) and *f?(0; w,y) can be completely expressed in terms of 
animal number generating function when wall is neutral i.e, in terms of A(l,y). We will 
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use these results in rest of the paper to study the surface effects for DBP in 2 and 3 
dimensions. 

In the adsorbed regime the number of monomers in direct contact with the wall is 
proportional to s and v a d = 0. This implies that the scaling form of 0(x, s) in the 
adsorbed regime would be 

^(x, a )~^(|x|/0 (18) 

where £ = (w — w c ) u is the chracterstic length scale in the system. Since we are away from 
the critical regime, £ is well behaved and never diverges for finite w. Also £ is independent 
of the size s of the polymers. The normalization of scaling function g(r) is chosen such 
that 

/oo 
cftq/flxl) = 1 (19) 
-oo 

A s (w)4>('K, s) is the coefficient of y s in the expansion of \l/(x; iu, y). In the adsorbed 
regime, A s (w) ~ (yoo(w))~ s for large s and behavior of 0(x, s) is given by Eq. (18), hence 
ty(x;w,y) will have a scaling form 

*(x;«;,y)~^(|x|/0 (20) 

where e = l-y/y^w). 

Since the scaling function g(|x|/£) has no y dependence, hence the scaling function of 
Cr(x; w, z) would also be just g(|x|/£) for w > w c . 

3 Two dimensional Directed Branched Polymer in presence of 
1-d penetrable surface 

For a penetrable surface, since the configurations spanning through the surface are al- 
lowed, there is no loss of entropy per monomer to take into account (Fig 1). Hence, 
w — 1 corresponds to a zero gain in free energy per monomer of the surface. We find that 
the value of y at which A(l,y) diverges, the only value of w which makes A(w,y) also 
divergent is w — 1. This implies that w c — 1 for a DBP in any dimension in presence of 
a Id line as long as A(l,y) gets singular at finite value of y. Then the polymer has bulk 
behavior at the critical point. At w — 1, i.e for directed branched polymer in bulk, we 
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have shown in an earlier paper [19] by scaling arguments and dimensional analysis 0(x, s) 
has a scaling form 



^(x^-s^/dx^) (21) 
This implies C (O, s) ~ s x ~ dVc , and the crossover exponent a is exactly given by 

a = \- dv c =\-Q (22) 

where v c is the transverse correlation exponent of a d+ 1 dimensional DA in bulk, which is 
equal to the correlation length exponent for a d dimensional HCLG with nearest neighbor 
exclusion. 

As we go to higher dimensions even though entropy loss and energy gain balances each 
other at w — 1, the polymer might start binding to a line only at wall activity greater 
than 1. For directed branched polymers, when A(l,y) has no divergence, w — 1 is not 
the critical point of the surface transition. Instead it is given by 

_ l + lM(l,y c ) 
Wc -l-y c /A(l,y c ) (23) 
where y c is the large polymer limit fugacity value of the polymer with neutral wall i.e, 
when w — 1. 

As an example, on a Bethe lattice with co-ordination number 3 the function A(l,y) is 

and y c = 1/4. At y — 1/4 the function As(l,y) = 2, and substituting in Eq.23 we get 
w c = 12/7, which is greater than 1. 

The 1 + 1 d DA gets mapped to a Id HCLG. For 1 + 1 dimensional DAs in bulk in 
[19] we have derived the exact expressions of A(l,y) and ty(x; l,y). Using them and Eq. 
14-16 we get the expressions for A(w,y) and ^(x;w,y) as follows 

A(w, y) = 2wyil + y) - (25) 

(1-y-wy- 3wy 2 ) + (1 + wy)^J (1 - 3y)(l + y) 

The connected density- density correlation function of the corresponding gas is simple 

exponential and hence the generating function ty(x;w,y) has a form 
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V(x;w,y) = K(w,y)exp(-b(y)\x\) 
where it is straightforward to calculate K(w,y) and b(y), and we get 



(26) 



K(w,y) 



2wy(l - 3y)(l + wy)(l-y + ^(1 - 3y)(l + y)) 
[(1 - 3y)(l + wy)y/T+jj +(l-y-wy- 3wy 2 )VT=3y] r < 



(27) 



and 



b(y) = log(/r+^ + y/l - 3y) - log(yi+^ - yfl - 3y) (28) 

The generating functions A(w,y) and ty(x;w,y) have branch cut at y — 1/3. For 
w — 1, they also have a pole singularity at y — 1/3. Hence, clearly the phase transition 
from desorbed to adsorbed phase occurs at w — 1, i.e w c — 1. For w < 1, y^w) = 1/3 
and for w > 1 it is given by the real positive solution of 



w - y - w(3 + 2w)y 2 - 3w 2 y 3 = (29) 
Near the critical point, for w = 1 + 5, to leading order we get yoo(w) to be 

1 5 2 

yoo{w) = 1- higher order terms (30) 

3 16 

The sticking fraction C st (w,y) can also be exactly calculated and we get it to be 



C st (w,y) 



y{l-w) 



+ 



1 + wy 



-i -l 



1 + 1/ ^(l + y )(l-3y)_ 
From this, near the critical point, we get the scaling form of C st (w,y) to be 



(31) 



C st {w,y) = y/l h{u) 
where e = 1 — y/y^w) and u = e~ 1 ^ 2 5 and we get 



(32) 



9m 2 
1+ 48" 



(33) 



This gives the order parameter C st (w) near the critical point to be proportional to y. 
For large values of w, expanding C st (w) in powers of (l/w) we get 
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Figure 2: Directed Branched Polymer on a Square lattice in presence of a Id impenetrable 
line about the diagonal. 



For w — > oo it approaches 1/2, the maximum possible fraction that can stick to wall, 
as expected. 

4 Two dimensional Directed Branched Polymer in presence of 
1-d impenetrable surface 

In the presence of an impenetrable surface, because of loss in entropy per monomer on 
the wall, the transition from desorbed to adsorbed phase takes place at a non trivial value 
of adsorption activity. 

Here we study a DBP in 1 + 1 dimension on a square lattice, in presence of an impen- 
etrable surface, about the diagonal (Fig 2). From the exact generating function A(l,y), 
A(w,y) and ty(x;w,y), it is straightforward to determine the critical value of w and the 
sticking fraction and the density profile in the desorbed, critical and adsorbed phase of 
the system. The 1 + 1 d case in presence of a solid wall along the growth direction can 
be mapped to the HCLG in 1-d with fugacity for all sites lieing along the negative axis. 
Making use of this mapping we get 




(34) 
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A(l , y)= (i^n*Em (35) 

and A(w,y) is easy to get by substituting in Eq. 14. This result can also be obtained 
using heap method. For an alternate treatment see [13]. 

The density-density correlation of the corresponding gas is exponential and hence the 
generating function \1/ (x; w , y) still has a form given by Eq. 26, but K(w,y) and b(y) are 
now given by 

K{ x = wyjl + wy)(y/l + y + y/l-'dy) 

{W> (l + y)(l-w 2 y 2 )^/T^^+(l-y-(A-w)wy 2 -w 2 y 3 )^/T+^ 1 ' 

and 

b(y) = log(vT+^ + yf\ - 3y) - M^TT^ - ^1 - %) (37) 

The generating functions A(w,y) and fy(x;w,y) have a branch cut at y = 1/3. At 
w — 1, has no divergence and y c = 1/3. Substituting in Eq. (23), we get w c = 3. 

This value is greater than the value for 1 + 1 d DBP with a penetrable surface. This is 
expected, since the tendency of polymer to grow away from the surface is more when the 
surface is impenetrable and hence only when the surface gets sufficiently attractive, the 
polymer starts sticking to it. For w > 3, the closest singularity to the origin occurs at 



y s = ^ (38) 

For w < 3 the branch cut singularity 1/3 dominates and hence y 00 (w), the infinite 
polymer limit fugacity value is equal to 1/3 for w < 3. Whereas for w > 3, yoo{w) = y s . 
Free energy is a constant and the order parameter, C st (w) is zero for w < 3. 

We get the sticking fraction, C st (w,y) to be 



C st (w,y) = 



l-2y-3y 2 + (-l + y + 2y 2 )J(l + y)(l- 3y) 



(39) 



y[-2y +{w + 2y- wy)^{\ +y)(l- 3y) + w(-l + 2y + 2y 2 )} 

Near the critical point for w = 3 + 5 and y = y s (l — e), we get the same scaling form 
for C st (w,y) as given by Eq. 32, with the scaling function h(u) to be 



1 u 
1 + 27 



(40) 
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Figure 3: Sticking Fraction in presence of line for a Directed branched polymer in 1+1 
and 2+1 dimensions, when the polymer size tends to infinity. 

Hence C s t(w) is proportional to y near the critical point and approaches 1/2 as 
w — > oo. This is plotted in Fig 3 along with the sticking fraction for the penetrable case. 
The qualitative behavior in both cases is just the same, the main difference being the shift 
of the transition point from 1 to 3 and the initial slope. For large value of w it is easy 
to expand C st {w) in powers of 1/w. It should be noted that large w expansion of C st (w) 
will involve powers of u> -1 / 2 in this case as well. 

Using the exact equations for generating function for ty(x;w,y) we translate these 
results to the constant number ensemble and we get the function (f>(x, s) in three re- 
gions giving the spread of sites as a function of distance from the wall. Here we give 
these calculations for the impenetrable case only because the qualitative behavior in both 
impenetrable and penetrable case is exactly same for 1 + 1 dimensional system. 

In the desorbed phase(u> < 3), expanding near y c = 1/3 as y = y c e~ e , we get the 
scaling form for ty(x; w, y) to be 

fy(x; w, e) = c(u>)exp(— xV3e) (41) 

where, c(w) is a w dependent constant and is 2(3- w) 2 • 

To obtain <j>(x, s) for large s, we need to determine the coefficient of y s in the series 
expansion of ^/(x;w,e) i.e, 

*(*; w,e) = Y^ 0(x, s)A s (w) y s = c(w) ^ \^-fr (1 - ^ (42) 
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= cW S (%) £ r[* + i] r[. + i]r[-*/2] (43) 

For fixed k and large s, 

r[s + i] v ; 

Hence the leading singular behavior of s)Ag(tu) in the desorbed phase is given by 

«*, .) A = — £ r[Jt + 1]r[ _ V2] (45) 

Since Y[—k/2] has poles when is an even integer, only odd terms contribute to the 
sum. It is easy to sum the resulting series, giving <f>(x, s) for large s in the desorbed phase 
to be 

<f) de {x, s) = exp (-^T^ ( 46 ) 

For w — 3, c(w) is singular and we have to keep terms till first order in e in the 
expansion(for w < 3 the constant term dominates) and we get 

3, e) = -exp(-xVSe) (47) 

Again, just as in the desorbed phase expanding fy(x; 3, e) in powers of y s , the average 
number of sites at a distance x, i.e <f>(x, s) for the critical region for large s is 

^, a ) = — erfc^j (48) 

Hence we see that at w = 3 not just the crossover exponent a is equal to 1/2, but 
even the scaling form of (f>(x, s) is same as that of a (1 + 1) dimensional DA in bulk [19] 
and hence same as that of the penetrable wall at the critical point. This unusual result 
can be understood as coming from exact cancellation of decrease in entropy and increase 
in internal energy at the critical point. Also note that the value of exponent a = 1/2 for 
DBP is equal to the estimates of a for branched polymers [21] and linear polymers [22] 
in 2 dimensions. Infact for adsorption of an undirected d dimensional branched polymer 
to a d — 1 dimensional surface, the crossover exponent a is conjectured to be 1/2 in all 
spatial dimensions [23]. 

14 



1 




Scaled transverse co-ordinate 



Figure 4: Denisty profile of a 2 dimensional DBP in presence of a ldimensional surface 

For w > 3, the behavior of the generating function is dominated by the singularity 
given by Eq. (38). For w ^> 3, y s ~ \j \fw and we get the large s behavior of (p(x, s) to 
be 

(f)(x, s) = s exp(— x) (49) 

i.e, most of the sites stick to the origin as expected 

Similarly, expanding A(w, y) about y c and then going to constant number (s) ensemble, 
we get A s (w) for large s as A s (w) ~ 2^pc(w)3 s s 2 in the desorbed regime. Hence the 
number of animals in presence of the Id impenetrable wall i.e A s (l) for large s are A s (l) ~ 
2^7=3 s s~2. This gives 9d e to be 3/2. This is consistent with the result derived for lattice 
trees by De'Bell et al [12]. Also we get at the critical point w — 3, A s (3) ~ "^^^ Ss_ ^' 
implying 9 C to be 1/2. For w ^> 3, A s (w) ~ (t^w) 13 , giving 9 ad = 0. 

The function <p(x, s) gives the density profile of the polymer as a function of distance 
from the surface. Since the configurations are very different in two phases as shown 
schematically in Fig 2 hence <p(x, s) is very different in three regions. In desorbed phase it 
peaks away from the surface at a distance of the order of the average transverse diameter 
of the polymer in the large s limit. Whereas at the critical point it peaks at the surface 
(Fig 4). 
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5 Three dimensional Directed Branched Polymer in presence of 
an attractive line 

In 2 + 1 dimensions, a DBP on a simple cubic lattice with nearest and next nearest 
neighbor connections gets mapped to the hard hexagon gas model in 2 dimensions at 
negative activity in the disordered regime, which was solved by Baxter [24]. He obtained 
the equation for the average density of the gas. It was shown by Joyce, there is an 
algebraic equation in z (activity of the gas) and p (density of the gas) [25]. 

The equation given by Joyce is quartic in z and 12th order in p. For convenience we 
will just rewrite it here [25] 

p{\ - p) 11 - (1 - p) 5 P 1 (p)z + p 2 (l - p) 2 P 2 (p)z 2 - p 5 P 1 (p)z' + p n (l - p)z 4 = (50) 
where 

P x (p) = (1 - 13p + 66p 2 - 165p 3 + 220p 4 - 165p 5 + 77 p 6 - 22p 7 ) 
P 2 (p) = (1 - I3p + 63p 2 - 125p 3 + 6p 4 + 401p 5 - 689p 6 + 476p 7 - 119p 8 ) 

The density p of the HCLG is just the negative of A(l,y) and z = — y/(l + y). It is 
straightforward to get an algebraic equation in A(l, y) as a function of y [26]. As A(w, y) 
is a simple rational function of A(l,y), y and w (see Eq.14), substituting p in terms of 
A(w,y), the grand partition function of the 2 + 1 dimensional DA in presence of an one 
dimensional line about the main diagonal of the lattice, we get a 12 th order polynomial 
equation in A(w,y), where the coefficients are functions of w and y. Explicit writing 
down the equation is rather tedious and is omitted. Since A(l,y) becomes singular for 
y = y c = 2/(9 + 5\^5), hence in presence of 1-dimensional line the polymer will undergo 
a desorption- adsorption transition at w — 1. For w < 1 the dominant singularity will 
be y c and yoo(w) = y c . For w > 1, at y^w), A(w,y) tends to infinity and at this point 
the coefficient of highest order term must be zero. Since we have a 12th order equation 
in A{w,y) hence by equating the coefficient of the 12 th order term to zero, we get a 
polynomial equation in y and w (Q(y,w) = 0) whose smallest positive real root would be 
y^w). This polynomial is 12th order in w. But we can find the root numerically. The 
free-energy is just log(y OQ (w)) and hence can be evaluated numerically. 
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In this case the expressions of A(w,y) and other generating functions are rather com- 
plicated and hence it is difficult to go to the constant size ensemble. But at critical point 
system behaves like bulk and since 9 = 5/6, by hyper-scaling arguments v c = 5/12 which 
implies that the crossover exponent a — 1/6 (Eq. 22). And by Eq. 10, the sticking 
fraction C st (w,y) ~ (1 — y/y 00 (w)Y J ^ 6 as w — > 1 + asymptotically. 

By solving Q(y,w) = we get yooiw) as a function of w. Near the critical point for 
w — 1 + 5, to leading order we get 

yoo(w) = y c (l — c5 6 + higher order term) (51) 

where c = 5(57) 5 with 7 = (13y/E - 25)/50. 

In the large polymer limit, for y very close to y c , A(l,y) has a scaling form 



A(l,y)=a (l-j- 



1 + 01 l 1 "^ J + 



(52) 



where a = (v 7 ^ 6 ) -1 [26]. 

Hence taking y = y 00 (w)(l — e) and w = 1 + 5 we get the scaling function of C st (w, y) = 
e 5 / 6 h(u) to be 

h(u) = — ^(l + oz 6 )§ -6cu 5 (53) 

where u = 5e~ l l % . The scaling function h(u) is a function of w and y, which are both 
thermodynamic variables. 

For large w, expanding in power of 1/tu we get 

~ - ib - li^l G;) (54) 

and 

1 1 / 3 Q 

C st H~---W- (55) 

As u> — > 00, C s t(w) approaches 1/2, the maximum possible fraction of adsorption. It 
is like order parameter of the surface transition. It is plotted in Fig 3. As is clear from 
the scaling function, the sticking fraction increases much slowly than in the 1 + 1 d case. 
This is expected as there the polymer in d + 1 dimensions was getting adsorbed at a d 
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dimensional surface whereas here a polymer in d + 1 dimensions is getting adsorbed on a 
d — 1 dimensional surface. 
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